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ABSTRACT 



o : 

f-j Aims. We spectroscopically characterize the Galactic Bulge to infer its star formation timescale, compared to the other Galactic components, 
CZJ ■ through the chemical signature on its individual stars. 

^ Methods. We derived iron and oxygen abundances for 50 K giants in four fields towards the Galactic bulge. High resolution (R=45,000) spectra 
J> | for the target stars were collected with FLAMES-UVES at the VLT. 

' ^ esu ^ s - Oxygen, as measured from the forbidden line at 6300 A, shows a well-defined trend with [Fe/H], with [O/Fe] higher in bulge stars 
r> [ than in thick disk ones, which were known to be more oxygen enhanced than thin disk stars. 

Conclusions. These results support a scenario in which the bulge formed before and more rapidly than the disk, and therefore the MW bulge 
can be regarded as a prototypical old spheroid, with a formation history similar to that of early-type (elliptical) galaxies. 

Key words. Galaxy: bulge - Stars: Abundances, Atmospheres 



1. Introduction 

The central regions of spiral galaxies, known as bulges, are 
made of stars in randomly oriented orbits for which two dis- 
tinct formation processes have been proposed. In the so-called 
"classical bulges", most stars originate in a short phase of star 
formation when the universe was only a few Gyr old. Instead, 
in the so-called "pseudobulges", stars form in the disk over a 
more extended period of time, and the bulge results from the 
secular evolution of the disk driven by the development of a 
bar (Kormendy & Kennicutt 2004). It is currently believed that 
the pseudobulge mode dominates in late-type spirals (Sc and 
later) and the old-starburst mode (similar to the formation of 
elliptical galaxies) dominates in early-type spirals (Sb and ear- 
lier). The Sbc Milky Way galaxy sits on the borderline, and 
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to establish the origin of its bulge (either old starburst or disk 
secular evolution), extensive spectroscopic observations have 
been undertaken to obtain precise element abundance ratios, in 
particular [Fe/H] and [O/Fe]. 

The detailed chemical composition of stars carries the sig- 
nature of the enrichment processes undergone by the interstel- 
lar medium up to the moment of their formation. Thus, elemen- 
tal ratios are sensitive to the previous history of star formation 
and can be used to infer whether there is a genetic link between 
different stellar groups. In particular, the relative abundances 
of iron and <?-elements play a key role because a-elements are 
predominantly produced by Type II supernovae (SNII) while 
supernovae of Type la (SNIa) dominate instead the production 
of iron. The SNIIs come from massive, short-lived stars, while 
SNIas result from binary evolution and are characterized by a 
very broad distribution of delay times (from a few 10 7 to over 
10 10 yr) between the star formation event and the SN explosion 
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depends on the timescale of star-formation and metal produc- 
tion (Matteucci & Greggio 1986). 

Among the a-elements, oxygen is especially important be- 
cause produced only by SNII's, and its forbidden line [OI] at 
6300 A is easily measured on the spectra of red giant stars. This 
line allows one to derive reliable abundances because its atomic 
parameters are well known; and since it is due to a transition 
from the ground state to a collisionally controlled higher level, 
it is not subject to non-LTE effects (Asplund et al. 2004). Based 
on this line, accurate measurements of the oxygen abundance 
in the Galactic thick and thin disk are also available (Bensby et 
al. 2004), while several efforts have been made to measure de- 
tailed abundances in the bulge (Mc William & Rich 1994, 2003; 
Rich & Origlia 2005). However, the bulge data were rather 
scanty, and the resulting evidence was somewhat contradictory. 
To overcome these limitations, high quality spectra for a large 
sample of bulge stars are presented. 

2. Observations and data analysis 



Table 1. Characteristics of the four bulge fields 



Nr. 


Identification 


/ 


b 


Rgc 


E(fi - V) 










(pc) 




1 


b = -6 Field 


0.21 


-6.02 


850 


0.48 


2 


Baade's Window 


1.14 


-4.18 


604 


0.55 


3 


Blanco b = -12 


0.00 


-12.0 


1663 


0.20 


4 


NGC 6553 Field 


5.25 


-3.02 


844 


0.70 



Spectra for a sample of ~ 1000 K giants in four bulge fields 
have been collected at the VLT-UT2 with the FLAMES fibre 
spectrograph. All the stars were observed with the GIRAFFE 
arm of the instrument, with resolution R~ 20, 000, while 58 of 
them have also been observed with the UVES arm, at higher 
resolution R~ 45, 000, in the range 5800-6800 A. A complete 
description of the whole GIRAFFE sample will be given else- 
where (Zoccali et al. 2006, in preparation), here we focus on 
the high-resolution UVES spectra, for which the S/N ~50 per 
resolution element has allowed us to obtain accurate oxygen 
abundances. The bluer half of the UVES spectrum has lower 
S/N and was only used to derive carbon abundances from a C2 
band at 5635 A. Observed targets include 1 1 stars in a low red- 
dening window at (I, b) = (0, -6), 21 stars in Baade's Window, 
5 stars in the Blanco field at (/, b) = (0,-12), and 13 stars in 
a field near the globular cluster NGC 6553 (Table [TJ In the 
colour-magnitude diagram, these stars are located on the red 
giant branch, about 1 magnitude above the red clump, with the 
exception of 13 stars in Baade's Window that instead are on the 
red clump itself. 

Individual spectra were reduced with the standard 
FLAMES-UVES pipeline, including bias, flatfield, interorder 
background subtraction, extraction, wavelength calibration, 
and order merging. All the spectra for each star (a number 
between 7 and 17, depending on the field) were then regis- 



gle spectrum per star. In each plate, one UVES fibre was allo- 
cated to an empty sky region. This spectrum was registered in 
flux to match the equivalent width of the sky emission lines of 
each object spectrum and then subtracted from it. The spectrum 
of a fast-rotating B star, thus containing only telluric absorp- 
tion lines, was coadded to itself using different radial velocity 
shifts, reproducing the shifts applied to the target spectra to re- 
move the heliocentric, variable component of the radial veloc- 
ity. Stars with telluric absorption blended with the O line were 
rejected. Stars with such absorption in the wings of the O line 
were included with a different symbol. The strong sky oxygen 
emission line prevents a measurement of the oxygen stellar line 
when the radial velocity of the star is close to zero, and these 
stars were eliminated from the sample. In the end, 8 stars were 
discarded, leaving us with 50 stars. 

All the stars have V,I magnitudes, obtained either from 
photometry performed on wide field WFI images from the ESO 
2.2m telescope at La Silla (fields 1,3,4) or from the OGLE cat- 
alogue (Udalski et al. 2002) (field 2). In addition, J, H, K s mag- 
nitudes are available for all of them from the 2MASS point 
source catalogue (Carpenter et al. 2001). The assumed redden- 
ing for each field is listed in Table while standard extinc- 
tion laws (Rieke & Lebofski 1985) were used to derive ex- 
tinctions in the other bands. The five de-reddened magnitudes 
were combined to obtain photometric temperatures from V — I, 
V - J, V - H, and V - K colours, according to the latest empir- 
ical calibration (Ramirez & Melendez 2005). The mean of the 
four values was used as a first guess for spectroscopic analysis. 
Photometric gravity was instead calculated from the classical 
relation 

T, M» 
logg* = log g + 4 log — + 0.4(M bol - M bolQ ) + log — 

1 Q Mq 

by adopting a mean distance of 8 kpc for the bulge, T =577O 
K, logg =4.44, M bo io=4.75, and M„=0.85 M . 

The equivalent widths for selected lines of Fe, Na, Mg, Al, 
Si, Ca, Sc, Ti, and Ni were measured using the new automatic 
code DAOSPEC (Stetson and Pancino in prep.). The selection 
of clean Fe lines, and their atomic parameters was compiled us- 
ing a spectrum of ft Leo as reference, observed at the Canada- 
France-Hawaii Telescope with the ESPaDOnS spectrograph, at 
resolution R=80,000 and S/N~ 500. The following parame- 
ters were determined for [i Leo: T e ff=4550 K, logg=2.3, mi- 
croturbulence velocity V t =1.3 km/s . Following Fulbright et al. 
(2006) we carried on our analysis relative to a nearby refer- 
ence star; i.e., we derived relative log gfs using Leo as a ref- 
erence, requiring [Fe/H]=+0.30 for each single line. More de- 
tails on the analysis of fj. Leo, the selection of clean lines, and 
the choice of its metallicity are given in Lecureur & al (2006). 
With the same atomic parameters, we obtain [Fe/H]=-0.52 for 
Arcturus and [Fe/H]=0.04 (e(Fe) G =7.55) for the Sun. This was 
then included as zero point for the iron determination of our tar- 
get stars. Atomic parameters for the [OI] and Nil lines around 
6300 A were instead assumed to be identical to those employed 
in previous measurements for disk stars, so as to allow a mean- 
ingful comparison between bulge and disk abundance patterns 
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Sun and // Leo, respectively. The damping constants were com- 
puted where possible, and in particular for most of the Fel lines, 
using the collisional broadening theory (Barklem, Anstee & 
O'Mara 1998, Barklem, Piskunov & O'Mara 2000). 

The LTE abundance analysis was performed using well 
tested procedures (Spite 1967) and the new MARCS models 
(Gustafsson et al. 2002). Spectrum synthesis was performed 
with turbospec (Alvarez & Plez 1998) and counterchecked with 
Barbuy et al. (2003), including the effects of molecular lines on 
the derived atomic abundances, which is of special importance 
in the case of oxygen, due to the formation of CO molecules, 
locking part of the oxygen, and TiO molecules, changing the 
continuum shape. Excitation equilibrium was imposed on the 
Fel lines in order to refine the photometric T e j, while photo- 
metric gravity was imposed even if ionization equilibrium was 
not fulfilled. The reason for this choice is that Fell lines are all 
very weak and contaminated by CN. On the other hand, while 
the strong and differential reddening has a strong effect on the 
photometric T e ff (why we prefer the excitation one), photomet- 
ric gravity is robust. In fact, an error of AE(B - V) = 0.05 
implies an uncertainty of A T e ff = IQQK in the photometric 
temperature but only Alogg=0.015 in the photometric gravity. 

A synthetic spectrum was calculated for each star, with the 
parameters and Fe given above, adopting solar abundance (see 
below) as a first guess and then iterating until the oxygen line 
was reproduced well. Other elements whose abundance might 
affect the oxygen were also measured in the iterations, namely: 
carbon from the C2(0,l) bandhead of the Swan A 3 n g -X 3 n„ 
system at 5635 A, nitrogen from the CN band at 6498.5 A, 
Ti, and Ni. The [C/Fe] ratio was found to range between -0.30 
and +0.10 for our stars. The accuracy we have in these mea- 
surements is not very high (+0.2 dex), but we verified that a 
variation of +0.2 dex in the adopted C implies a negligible 
variation in [O/Fe], except for the coolest stars at super-solar 
[Fe/H], for which it can introduce up to a 0. 1 increase in [O/Fe]. 
Nitrogen, on the other hand, ranges between [N/Fe]=0.00 and 
[N/Fe]=+0.50. Particular care was taken in measuring Nickel, 
given that the [OI] 6300.304 A line is blended with the two iso- 
topic components of the Nil line at 6300.335 A (Johansson et 
al. 2003). Thus isolated Ni lines were used to measure [Ni/Fe], 
and this value was assumed in the deblending of the [OI]6300 
A line. Standard solar abundances (Grevesse & Sauval 1998) 
were adopted as reference zero point for the abundances of 
bulge stars, except for the value of solar oxygen, where e(0) = 
8.77 was assumed (Allende Prieto, Lambert & Asplund 2001). 

From the dispersion of Fe lines, we estimate a statistical 
error of 0.05 dex on Fe. On the other hand, a systematic error of 
- 100K in T eff implies A[Fe/H] = -0.04 and a negligible effect 
on [O/H], while A[Fe/H] = +0.04 and a A[Fe/H] = +0.1 
would be given by a +0.3 dex increase in log g. 

A table with the complete list of parameters and abun- 
dances, including sodium, magnesium, and aluminum for each 
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Fig. 1. Upper panel: The [O/Fe] vs. [Fe/H for our bulge stars 
(circles) along with previous determinations in other bulge stars 
from optical (open triangles) and near-IR spectra (crosses). 
Open symbols in our measurements refer to spectra with lower 
S/N or with the O line partially blended with telluric absorp- 
tion. Lower panel: oxygen/iron trend in our bulge stars vs. that 
for thick and thin disk stars. The solid line shows a linear fit to 
the thick disk data points with [Fe/H]> -0.5 and is meant to 
emphasize that all bulge stars with -0.4 <[Fe/H]< +0.1 are 
more oxygen-enhanced than thick disk stars. This plot enlight- 
ens the systematic, genetic difference between bulge and disk 
stars, thus excluding that bulge stars were once disk stars that 
then migrated inward to build up the bulge. 



3. [O/Fe] ratios 

The resulting [O/Fe] vs. [Fe/H] plot is shown in Fig. ^ The 
upper panel shows the comparison between our measurements 
and recent oxygen abundance determinations for bulge stars. 
Earlier measurements (McWilliam & Rich 1994) are omitted 
here as more uncertain due to their much lower spectral res- 
olution. A few corrections have then been applied to ensure 
consistency between our bulge abundances and those from pre- 
vious determinations for bulge and disk stars. For the bulge op- 
tical data (McWilliam & Rich 2003), a shift of -0.06 in [O/Fe] 
takes the difference in the adopted solar oxygen abundance into 
account, their e(O) =8.71 versus our 8.77, while a further +0. 1 
in [O/Fe] corrects for the different \oggf--9.1\l for the [OI] 
line, 0.1 dex lower than the one adopted here. Similarly, the 
bulge abundances from near-IR measurements (Rich & Origlia 
2005) were shifted upwards by +0.06, due to their assumption 
of a higher solar oxygen, e(O) =8.83. This figure shows that 
our measurements do not conflict with previous results, but the 
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The lower panel of Fig. ^ shows the main result of the 
present investigation: the [O/Fe] vs. [Fe/H] ratios in the bulge 
compared with those for the thick and the thin disks (Bensby 
et al. 2004). These measurements are as consistent as possi- 
ble, in the sense that they come from the same line, and we 
have intentionally adopted the same atomic parameters both for 
oxygen and for nickel. The solar oxygen abundance to which 
the disk stars were referred was lower by 0.06 dex, so an iden- 
tical downward shift was then applied to their measurements. 
This plot shows that the thin disk, thick disk, and bulge evolved 
through different chemical trajectories. In other words, bulge 
stars did not originate in the disk and then migrate inward to 
build up the bulge, but rather formed independently of the disk 
(Minniti 1995, Ortolani et al. 1995). Moreover, the chemical 
enrichment of the bulge, hence its formation timescale, has 
been faster than that of the thick disk, which in turn was faster 
than that of the thin disk (Matteucci, Romano & Molaro 1999). 

A few stars in our sample show lower [O/Fe] than the mean 
locus of thick disk stars. Their number is consistent with the 
expected number of foreground disk contaminants. According 
to the Besangn model (Robin et al. 2003) and to the disk con- 
trol field in Zoccali et al. (2003) the contamination fraction is 
15% in all our fields except the Blanco one, where it is 45%. 
All except two of these stars would indeed require stronger 
gravity according to the ionization equilibrium. Two extreme 
cases (logg~3-4) have been excluded from the sample so are 
not shown here. We also excluded stars with radial velocity 
close to zero (with a sky emission line on top of the oxygen 
line), further biasing our sample against disk contaminants. We 
therefore conclude that the stars above the thick disk line in 
Fig. 1 are likely to be bona fide bulge stars. 

It is evident from Fig. 1 that the dispersion in [O/Fe] for 
bulge stars is higher than for both thin and thick disks. While 
at least part of this effect is expected, given that disk stars are 
nearby, bright and isolated stars for which higher S/N spectra 
were obtained, the possibility that part of the dispersion is real 
cannot be excluded, given that in the bulge we are sampling 
stars in a region much larger than the solar neighbourhood. It is 
worth mentioning that no [O/Fe] gradient was present among 
the four different fields. 

4. Conclusions 

Age-dating of the bulge stars from their colour-magnitude dia- 
gram has already demonstrated that the bulk of them are older 
than ~ 10 Gyr, with no detectable trace of intermediate age 
stars (Ortolani et al. 1995; Kuijken & Rich 2002; Zoccali et 
al. 2003). Interpretation of the evidence now provided on the 
[O/Fe] ratios is that the formation of the Galactic bulge was 
faster than that of the disk and preceeded it. The most metal- 
rich bulge stars were the latest formed during the bulge evolu- 
tion, and their ages correspond to the exhaustion of interstellar 
matter in the bulge. The present results reveal that the oxygen- 
to-iron starts to drop at the metallicity of [Fe/H]«;-0.2 due to 
the progressive enrichment by SNIa. A fast star formation rate 
is needed to explain this result. 
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Fig. 2. Comparison between the observed [O/Fe] ratio of bulge 
stars and the prediction of the simulations by Immeli et al. 
(2004). See text for details. 

borhood. Strictly speaking, our result implies that the bulge 
stars did not originate in disk material (either already stellar 
or gaseous) coming from as far as the solar neighbourhood. 
Due to the possible presence of radial gradients in the disk, a 
comparison with stars in the inner disk would have been more 
conclusive. Unfortunately no such sample is available at the 
moment. An alternative approach to interpreting of these re- 
sults is the comparison with theoretical simulations. In Fig. [2] 
we compare our results with the predictions of bulge [O/Fe] 
ratios by Immeli et al. (2004). They simulate the formation of 
galaxies from clouds with different dissipation efficiencies. The 
case with high dissipation (model B) rapidly forms massive 
clumps that spiral to the centre and merge to a central bulge 
component in a strong starburst. In contrast, when the dissi- 
pation efficiency is low (model F) the galaxy forms a thicker 
disk, and only at later times does instability set in, driving the 
formation of a bar and then a pseudobulge. Model D is the in- 
termediate case. Clearly, our data can exclude model F, while 
it agrees with model B very well. This makes the MW bulge 
similar to early-type galaxies, in being a-element enhanced, 
dominated by old stellar populations, and having formed on a 
timescale shorter than ~ 1 Gyr (Thomas et al. 2005). Therefore, 
like early-type galaxies the MW bulge is likely to have formed 
through a short series of starbursts triggered by the coalescence 
of gas-rich mergers, when the universe was only a few Gyr old. 
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